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Abstract 

Stochastic representation for interaction of quan- 
tum systems is formulated which allows to re- 
place some of them by equivalent but purely 
commutative random sources. The formalism 
is applied to two-level systems interacting with 
Gaussian thermal bath. Strong-coupling non- 
Marcovian effects and besides long-living fluctu- 
ations in common susceptibility of two systems 
subjected to the same bath are considered. 

1. Stochastic representation. 

In quantum statistical physics and kinetics 
typically one deals with tasks about interac- 
tion between some "dynamical system" (DS, mi- 
croscopical or with only a few degrees of free- 
dom) and its macroscopic surroundings ("ther- 
mal bath", or "thermostat"). Usually, an inter- 
action Hamiltonian has bilinear structure (or can 
be reduced to that): 



Hint = X! * ' 



(1) 



where operators Dj and Bj act in different lin- 
ear (Hilbert) spaces D and B (native spaces 
of DS and thermostat, respectively). In the 
Heisenberg picture, the operators Bj(t) = 
exp(iHht)Bj exp{—iHht) (with Hf, being Hamil- 
tonian of autonomous thermostat) play like ran- 
dom perturbations with respect to DS. But, 
in opposite to classical Langevin forces, they 
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are operator-valued non-commutative random 
sources and therefore withdraw the DS evolu- 
tion into the immense direct product of the two 
spaces, D ® B. 

Possibly, the task would become more sim- 
ple, if Bj{t) were replaced by equivalent commu- 
tative random sources, so that the DS evolution 
would remain in D . The equivalence means ex- 
act reproduction of results under interest, first 
of all, exact simulation of role of the thermostat 
temperature and and effects of DS self-action 
mediated by thermostat, in particular, dephas- 
ing and dissipation. 

The variant of such a replacement was sug- 
gested in [1,5]. Let Hd^t) be Hamiltonian of free 
DS. Introduce the notation 

H{B,t) = Hdit) + Y.B,it)Dj (2) 

j 

This is joint Hamiltonian in the representa- 
tion which is Heisenbergian with respect to 
the thermostat. Let some DS observables, Jk, 
are staying under actual observation (measure- 
ment). Then statistical operator of complete sys- 
tem "DS plus thermostat", R{t) , undergoes the 
equation 

R = v{t)J oR + i {RH{B, t) - H{B, t)R} , (3) 

where v(t)J = J2''^k{t)Jk , with Vk{t) being test 
functions of the observation and o the symbol 
of symmetric product (Jordan product). At 
Vk (^) = , this is the usual von Neumann equa- 
tion for joint density matrix. At Vk{t) ^ , the 
trace of R{t) m. D ® B represents [1,5] the char- 
acteristic functional of the observables, H(i,i;): 



- [ v{t')J{t')dA X U{t, /) = -i{Hb + Y: mBj}U{t, /) 

2 J J .7 



xfecp 



(4) 



Here Rin = R{—oo), the presence of time ar- 
gument in Jk{t) means that these operators 
are treated in the Heisenberg picture, left-hand 
(right-hand) oriented arrow indicates chronolog- 
ical (anti-chronological) time ordering, and the 
integrals are taken over the region t' < t (with 
t being present time moment). Let us replace 
the operators B{t) in Hamiltonian (2) by com- 
mutative (similar to C-numbers) variables ^(t) = 
x{t)+iy{t)/2 or ri{t) = x{t)—iy{t)/2 and, instead 
of Eq.3, consider the equation 



Clearly, the time-dependent operators Bj{t,f) 
characterize non-autonomous behaviour of the 
thermostat. On the left side of Eq.7, g{t) and 
f{t) are test functions for the random sources 
x{t) and y(t), respectively. On the right side, 
according to Eqs.7 and 8, g{t) arc test functions 
for the thermostat obscrvablcs B{t) while fj{t) 
are classical (C-number valued) forces which per- 
turb the thermostat being conjugated with Bj{t) 
(in the sense of non-equilibrium thermodynam- 
ics). Note that the relation (7) can be rewritten 
in the form [1]: 



p{t) = v{t)J o p{t) + L{t)p{t) ^ 



(5) 



exp J J^MMt) + fj{t)yjit)]dt\ = 



L{t)p = i{pH{rj,t)-H{C,t)p}= (6) = TVs e5q5 J 1 5^ 



9 jit) 



+ if jit) 



Bjit)dt 



= Y,yj{t)Dj o p + i[p,H{x,t)] , 
j 

treating all the 77 (t), ^(t), x{t) and y{t) as ran- 
dom processes which reflect an influence by the 
thermostat. Assume that Rin has the factorized 
form: Rin = p{—oo)ph, where p5 = TrD^( — 00) 
is statistical operator of the thermostat. All of 
x{t) and y(t) are defined if their characteristic 
functional is defined. We make it as follows: 

^exp J [g{t)xit) + f{t)y{t)]dt^ = (7) 



xfecp 



E 



9 jit) 



-ifjit) 



Bjit)dt\ pb. 



demonstrating formally equal rights of x(i) and 
yit) on the right-hand side of Eq.7. 

As it was shown in [1], if statistics of com- 
mutative random sources in Eqs.5 and 6 is deter- 
mined by Eqs.7 and 8 then the sources exactly 
imitate real quantum thermostat. Concretely, 



(p(t)) =IVsi?(t), 



(9) 



= TVBe5q^|^y git)BitJ)dt^ x 

x^{\J 9it)Bit,f)dt^ Pb 

(the subscripts are omitted for brevity), where 
the operators Bjit,f) are defined by means of 
the equations 

BjitJ) = U+itJ)BjUit,f), (8) 



Ttd {pit)) = Eit,v), (10) 

where the angle brackets (...) denote statistical 
averaging with respect to x(t) and yit) (or to 
rjit) and ^(t)) in accordance with (7-8). 

2. Properties of the random sources. 

The pay for the desired replacement is the 
doubling of random sources: instead of each real 
observable (Hermitian operator) Bjit) wc obtain 
the pair of either real variables (xj(t) and yjit)) 
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or mutually conjugated complex variables (rjj 
and = rjj)- From Eqs.5-6 it is seen that Xj(t) 
play direct role of random forces (potentials), 
while yj{t) stay in positions of test functions 
which correspond to observing DS by thermostat 
(in analogy with Vk(t) in Eq.3). In absence of 
yj{t), under random pump introduced by Xj{t), 
both the energy and entropy of DS would grow 
as much as possible. However, like in any mea- 
surement process, influence of yj{t) introducing 
observation by thermostat to the "stochastic Li- 
ouville operator" (6) destroys unitarity of evolu- 
tion of the statistical operator p{t) and decreases 
phase volume (entropy) of DS. At the same time 
yj {t) are responsible for the energy outflow back 
to thermostat, i.e. for dissipation, and hence 
for non-uniform (thermal) DS energy distribu- 
tion over states of DS. It should be emphasized 
that according to (9-10) on the average the uni- 
tarity is ensured. 
Due to Eqs.7-8, 




Here from it is clear that {y{t)) = 0, {y{t)y{t')) 
= 0, and all the higher self-correlators of y{t) also 
are equal to zero. But the cross-correlators be- 
tween y(t) and x{t) differ from zero representing 
the response of thermostat to its perturbation 
by DS. Thus, y{t) are not C-numbers in literal 
sense. But the singularity of statistical proper- 
ties of y(t) only facilitates factual calculations. 
Importantly, y{t) can possess correlations with 
more late values x{t' > t) only, in correspon- 
dence with the causality principle. 

In case of Gaussian equilibrium ("bosonic") 
thermostat, all the functions (11) can be ex- 
pressed through pair correlators 

f°° duj 

{Xj{T)XmiO)) = COs{uJT)Sjm{uj) , (12) 

Jo TT 



{xj{T)ym{0)) = 2^{t) / sin(a;T)x 
Jo 

, / a; \ ^ , -.duo 
xtanhl — 1 Sjm{uj) — 

Here T is temperature of thermostat, i?(T) 
is Heavyside function, and Sjm{<^) is a non- 
negatively defined spectral matrix. The formu- 
las (12) show that neglecting the y(t)-component 
of random sources in (5-6) in essence would be 
equivalent to infinity of thermostat temperature. 

In Eqs.9-10, only the average value of the 
"stochastic density matrix of DS" , p{t) ,does ap- 
pear. Naturally, its higher statistical moments 

{p{t)^...®p{t)) (13) 

describe several copies of DS interacting with the 
same thermostat. The undistinguishness of the 
copies require to concretize their quantum statis- 
tics (second quantization rule). In [1] we con- 
sidered how the case of Fermi statistics reduces 
to analysing the moments (13), thus replacing 
second quantization in a "many-particle prob- 
lem" by equivalent "second randomization" in 
one-particle problem. In this approach, an addi- 
tional direct (e.g. Coulombian) interaction be- 
tween samples of DS can be included as interac- 
tion through the second (Gaussian) thermostat, 
i.e. additional pairs of stochastic sources x{t), 

y{t). 

3. Two-level system in Gaussian ther- 
mostat. 

Let us apply the above formalism to two-level 
system (TLS) interacting with Gaussian thermal 
bath. We can put on 

and H(i{t) = u{t)J/2 . Here u{t) is energy differ- 
ence between two states, possibly depending on 
time because of external perturbation of TLS. 
The operator D performs contact with ther- 
mostat which forces random switchings of TLS 
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state. The operator J corresponds to obser- 
vation of TLS state, for instance, spin orien- 
tation, or alternate velocity of quantum parti- 
cle in the simplest model of Brownian motion. 
We use this model to illustrate role of the y{t)- 
component of random source, i.e. "dissipative" 
term y{t)D o p{t) in Eqs.5-6. 

First, let both states of free TLS have equal 
energies: u{t) = 0. Then the Eqs.5-6 can be 
integrated up to rather visual expressions. In 
particular, for the trace 

e{t, v) = Tr£) p{t) = pii + P22 

one can derive the recurrent equation 

e{t,v)=Y{t,-oo)+ (15) 

+ / Y{tM)v{ti)X{tiM)v{t2)e{t2,v)dtidt2, 



X{ti,t2) = COS 



Y{ti,t2) = cosh 



t2 



x{t')dt' 



t2 



y{t')dt' 



where t> ti> t2, and E{t,v) = {9{t,v)). 

Consider the correlation function of spin fluc- 
tuations (or velocity fluctuations if speak about 
Brownian particle, or so on), K{t) = (J(t)J(O)), 
and the diffusivity A = /q°° K{T)dT. First of all 
the case is curious when the noise of thermostat 
is "white", i.e. S{ijj) = const in (12). To be pre- 
cise, it should be underlined that both the cor- 
relators (12) can not simultaneously turn into 5- 
functions, and in this sense quantum noise never 
can be white. After expanding Eq.l5 into series, 
averaging it with taking into account Eqs.l2 and 
then differentiating by the test function, one ob- 
tains 



K(t) 



e"'^^'^ cos 



2S_ 



Intanh I — - — ] dr 
V 2 



Under increase of the order of magnitude of 
the dimensionless noise intensity, S/T, initially 
monotonous relaxation changes to more and 
more oscillating relaxation. Prom the formal 
point of view, this is just the effect of the y{t)- 
component, i.e. inverse influence by DS onto 
thermostat (indeed, in absence of y{t) the cor- 
relation would be purely exponential indepen- 
dently on S/T). 

Howevere, practically the opposite approxi- 
mation can occur more appropriate: 

S{u:) = S/[l + {Touf] , roT » 1 , (17) 

when sufficiently large correlation time of the 
thermostat noise (characteristic time of ther- 
mostat response), tq, excludes high frequencies. 
Then it is convenient to introduce the coupling 
energy e by means of 5(0) = 2e^ro . Calcu- 
lating the correlation function under assump- 
tions (17), we find again that at weak coupling 
(e/T < 1) monotonous (nearly exponential) re- 
laxation takes place, while at e/T > 1 (strong 
coupling) 



K{t) k exp(-2eV^) cos(eV/r) , 



(18) 



(16) 



that is again the oscillations arise and then be- 
come multiplied. 

Clearly, oscillations of relaxation lead to an 
excess decrease of the diffusivity A (spectral 
power density of J{t) at zero frequency) as com- 
pared with its value Aq which would realize in 
absence of y{t). In the case of Lorenz noise with 
Ttq > 1, the diffusivity can become extremely 
small resulting in spatial localization of Brown- 
ian particle. 

The non-monotonous relaxation means that 
the TLS response to periodic external pertur- 
bation can acquire "resonant" character having 
maximum at some non-zero frequency of the per- 
turbation (or, if the frequency is fixed, maxi- 
mimi at some non-zero level of thermostat noise). 
Such a kind of phenomena is known as "stochas- 
tic resonance" (see, for example, the review [2]). 
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We would like to emphasize also that oscil- 
latory relaxation represents non-Marcovian ef- 
fect which could not be adequately considered 
by the theory of Marcovian quantum dynami- 
cal semi-groups, i.e. time-local kinetic equations 
(this theory takes beginning from the classical 
work [3]). In terms of this theory, we considered 
"stochastic dilation" of a semi-group which is as 
much non-Marcovian as strongly the correlator 
(x(T)y(O)) determining oscillations in (16) and 
(18) differs from 5-function. 

4. Stochastic linear response of TLS to 
non-stochastic perturbation. 

Because of the symmetry (degeneration) of 
states of TLS, the trace (15) contains even pow- 
ers of v{t) only. If u{t) ^ then this symmerty 
breaks down. Define integral operators /, C and 
S by the formulas 

if{t)= f f{t')dt', 

J —oo 

C, S f{t) = cos, sin I j^* u{t")dt"^ f{t')dt' , 

At u{t) ^ instead of (15) the following equation 
can be derived: 

e = l + iyCye+ (19) 

+i{v + 2ySx)[l + dxCx\-^i{v + 2xSy)e 

Now, in the expansion of 9 into power series of 
v{t), 

e{t,v) = eo{t) + f j{to)v{to)dto+ 

J —oo 

+ \ f f K{h,t2)v{ti)v{t2)dtidt2 + ... , 
^ J —oo J —oo 

the first order term is present (as well as other 
odd terms). For simplicity, consider this contri- 
bution in the "hot thermostat limit": S/T « 
1 , Ttq >> 1 , when in accordance with (12) any 



multiplier y{t) contributes excess order of small- 
ness. Besides, assume the perturbation u{t) is 
infinitely weak. Then, saving lowest powers of 
y{t) and u{t) only, we can write 

J{to) = J+{to) + J-{to), (20) 

J+ = j^" sin|2^*°a;(i')di'| X 
xn(ti) / y{t2)dt2dti , 

J —oo 

J-= f sin (2 / ' x{t')df'\u{ti) f y{t2)dt2dti 

Of course, here the factual observation time to is 
more early than formal running time t. Under 
this condition, t drops out from any result of sta- 
tistical averaging, therefore it may be convenient 
to set t = 00. 

Obviously, J{to) disappears if we neglect 
y(t)-component of the bath noise. Consequently, 
this component (in company with x{t)) is re- 
sponsible for non-zero response of DS to exter- 
nal perturbations. The integral operators which 
act onto u[t) in the expressions J+ and J_ are 
nothing but random linear response functions 
(random susceptibility, or mobility, or so on), 
in other words, these operators represent mul- 
tiplicative (parametric) noise. What is hidden 
beyond it is the randomness of quantum proba- 
bilities of quantum transitions, i.e. dependence 
of the probabilities on thermostat noise [4] . 

Formulas (20) clearly highlight that the re- 
sponse functions in their non-averaged form has 
no decay and can connect past and fiiture events 
regardless of how far they arc distanced in time. 

5. Two TLS in shared thermostat. 

Consider two identical TLS contacting with 
one and the same thermostat. Suppose that the 
systems are formally distinguishable and that 

they give additive contributions to the variable 
under observation. Then characteristic func- 
tional of the summary observable is presented 

by 
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(21) 

instead of (10). For its correlation function, to 
be designated as K2, in the hot thermostat hmit 
the Eq.21 yields 

K2{h,t2) = 2{(K{h,t2)) + (J+(tl)J_(i2)) + 

+ {Mh)J+{t2)) - 2 (J+(ii)) {J+{t2))} , (22) 

where ti > t2 is mentioned. 

Wc do not discuss the first term on the right- 
hand side of Eq.22. What is for the next terms, 
they describe "excess noise" of summary observ- 
able which is proportional to squared perturba- 
tion and comes just from the fluctuations in 
random linear responses (similar in both TLS 
and determined by Eqs.20). If ti — ^2 >> '^o 
where Tc is correlation time of the fast com- 
ponent of the J{t) fluctuations (that is decay 
time of equilibrium correlation function K(t) at 
u(t) = 0), then a contribution from the cross- 
correlator of J+ and J_ does survive only, and 
further it does not decay at all. A rough for- 
mal estimate of these long-living fluctuations by 
order of magnitude gives 

|(J+(oo)J_(-oo))|~(T2/ro^)(j)' 

Here the ratio Tc/tq crucially depends on the pa- 
rameter ero and can be much greater than unit 
(at ero << 1) as well as much smaller than unit 
(at ero >> 1)- 

Interestingly, the expression J- (to) in (20) 
corresponds to "reaction of the past on the fu- 
ture" because it includes future values of all the 
variables only with respect to the observation 
time. With no doubts, any correlator in which 
time argument of J_ is most late should turn into 
zero, otherwise it would have no physical sense. 
However, non-zero correlations between J_ and 
more late J+ does not contradict to causality 
principle. 

The non-decaying tails of such correlations 
in two-TLS system are quite similar to proper- 
ties of fluctuations of quantum transport prob- 
abilities in many-electron system considered in 



[4]. Their physical origin is the coherence (uni- 
tarity and time reversibility) of joint evolution 
of DS and thermostat: the coherence can be in- 
terrupted by external observation acts but con- 
serves between them. What is important, un- 
like [4], in the present example formal analysis 
of non-decaying correlations is not restricted by 
any time frames. 



6. Appendix. Generalized stochastic 
representation. 

The text above almost coinsides with the pa- 
per in Russian [5] (to be translated into English 

in JETP Letters). Since the proof of the ba- 
sic relations (9), (10) is omitted in [5] (see [1] for 
it), below we compensate this by proving gen- 
eralized stochastic representation which includes 
the above one as a particular case. 

Let the joint evolution (Liouville) operator, 
L , of combined system "DS plus thermostat" has 
the bilinear form as follows: 



L = Ld + Lb + ^ AojA^j 



(23) 



where Ljj and Lb are evolution (Liouville) op- 
erators of autonomous DS and thermostat, re- 
spectively, and Adj and Abj are some superop- 
erators acting in D and B , respectively (to be 
precise, they act in the spaces of linear operators 
defined in D and B). For instance, the super- 
operators Aoj and Abj can be Liouville opera- 
tors (commutators) or multiplication operators 
(Jordan products). For shorteness, assume that 
the terms X^ffc(i)Jfc corresponding to external 
watching for DS are included into Lo- 

Notice, that, from one hand, if a{t) is any 
random process with characteristic functional 



{exp 



f{t)a{t)dt 



(24) 



(angle brackets mean averaging over a's proba- 
bility distribution) , then the average of any func- 
tional ${«} can be formally represented by 

8 

da 



(25) 
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Similar relation holds also for a set of random 
processes aj(t). 

From another hand, since all oi Lq and A.Dj 
commute with all oi Lb and A^j, we can write: 



Bj 



dt 



(26) 



6xp 



Bj 



6aj{t) 



dt 



dt} X 



a{t)=0 



Here mutual chronological ordering of the two 
chronological exponents in their product is taken 
in mind. This is merely the consequence from 
chronological ordering and besides from the for- 
mal identity 

exp(0i02) = exp (oi-^-^ exp(a02) |a=o , 

where Oi and O2 are arbitrary mutually com- 
muting objects: [Oi, O2] = . 

Suppose that initial statistical operator i?j„ 
factorizcs, = p(— 00) pf, where p stands for 
DS density matrix. Then, comparing formulas 
(25) and (26), we conclude, firstly, that from the 
point of view of DS its joint evolution together 
with thermostat is equivalent to its individual 
evolution but influenced by the set of random 
sources aj{t) and governed by stochastic evolu- 
tion operator as follows: 

p = {v{t)Jop + LD + J2^j(*)^Dj}p (27) 

(we extracted external DS observation back from 
Ld). Secondly, in analogy with (24) and (25), 
characteristical functional of these sources is de- 
termined by 



Formulas (27) and (28) give the stochastic 
representation of the joint evolution. In fact, 
it is obtained by nothing but a variant of the 
Stratonovich transformation. Dependently on 
concrete contents of operators Adj and Agj the 
random sources aj{t) can behave either simi- 
lar to usual classical random processes (like x{t) 
above) or in a singular manner (like y{t) above) 
or in some mixed fashion. Of course, in general 
all the sources are mutually correlated. 

In particular, the above considered case of 
the bilinear interaction Hamiltonian (1) can be 
reduced to the representation (27)- (28) or de- 
rived from it if take into account the identity 

[R, DB] = [R,D]oB + [R, B]oD (29) 

which is valid for any pair of mutually commut- 
ing D and B , at [D,B] = . Hence, every pair 
Dj and Bj from (1) produces two terms in (23). 
Clearly, in the corresponding pair of a's one con- 
jugated with the first term on the right-hand side 
of (29) plays the role of x{t) (see Sec.l), while an- 
other (conjugated with second term) plays the 
role of y{t). 

It is useful to demonstrate a simple case 
which does not reduce to bilinear interaction (1). 
Concretely, let the joint evolution operator L is 
defined by 

LR = {Ld + Lb)R + iD' o [R, B] + iB' o [R, D] , 

(30) 

where D' differs from D, B' differs from B, and 
Ld,bR = i[R-,H(i^i,]. Then, instead of the repre- 
sentation expressed by Eqs.5-8, we arrive to the 
representation 

p = v{t)J op + y{t)D' op + i[p, Hd + x{t)D] , 

(31) 



exp 



YMt)a,{t)dt 



(28) 



exp / [git)xit) + f{t)y{t)\dt 



(32) 



TVs gxp |y [Lb + ^ fj{t)kBj\ dtj 



Pb 



= Trs exg 



g{t)B'{tJ)dt^ 
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xSxpj^y g{t)B'it,f)dt^ pb , 

B'{t,f) = U+{t,f)B'U{t,f) , (33) 

UitJ) = -i{Hb + f{t)B}U{tJ) 

Now the two random sources involve not two 
but four operators D, D', B and B' . Neverthe- 
less, again it is possible to separate analysis of 
thermostat, under given classical perturbation, 
and analysis of DS under influence by commuta- 
tive stochastic sources with given statistics (and 
again statistical properties of a number N of DS 
copies interacting with the same thermostat can 
be obtained from the moments < (Tr^ p)^ >, if 
the copies are distinguishable and undergo Boltz- 
mann statistics, or as it was considered in [1], if 
the copies undergo Fermian second quantization 
statistics). Concrete application of this formal 
example will be done elsewhere. 
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